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Abstract 

The electromagnetic form factor of the r\ meson is studied using 
the Light Cone QCD sum rules approach taking in to account higher 
twists contributions and SU{3) flavor symmetry breaking effects. 
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1 Introduction 



The QCD sum rules is one of the most powerful tools in studying the 
low energy hadron physics on the basis of QCD. In this method, physical 
observables of hadrons are related with QCD vacuum via a few condensates. 
In the literature, a new, widely discussed alternative to traditional sum rules, 
namely light cone QCD sum rules (LCQSR) is a convenient tool for the study 
of exclusive processes which involve the emission of a light particle. 

The LCQSR is based on the operator product expansion on the light 
cone, which is an expansion over the twist of the operators rather then the 
dimensions as the traditional QCD sum rules. In this method all the non 
perturbative dynamics encoded in the so called wave functions, which repre- 
sents as the matrix element of the nonlocal operators between vacuum and 
the corresponding hadron (more about this method can be found in pHIHl)- 
One of the promising ways for obtaining information about wave functions 
is experimental study of form factors. The pion electromagnetic form factor 
was calculated in the framework of LCQSR in jHE] • Recently the electromag- 
netic form factors of pion and Kaon is analyzed in [6] taking into account the 
effects of higher twist contributions and the SU(3) flavor symmetry breaking. 

The aim of this work is the calculation of the electromagnetic form factor 
of the rj meson taking into account twist 3 and twist 4 wave functions of the 
rj meson as well as the SU(3) symmetry breaking effects, i.e. the mass of the 
strange quark and the different value of the strange quark condensate. 

The paper is organized as follows: In Sect. 2, we derive the sum rules for 
the electromagnetic form factor of the r\ meson which is then analyzed nu- 
merically in section 3. Summary of our results and conclusions are presented 
at the end of Sect. 3. 

2 Sum Rules for the r\ Meson Electromag- 
netic Form Factor 

In this section we calculate the e.m. form factor of the r\ meson employing 
LCQSR. For this purpose, we consider the following correlation function: 
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where 

3u = 77g (ujulsu + dj u j 5 d - 2s7„7 5 s) (2) 
is the 77 meson interpolating current and 

it= E ( 3 ) 

is the electromagnetic current. 

Here, we would like to remark that in this work we will follow the usual 
parametrization of the 77 — rf mixing in octet-singlet basis and due to the 
smallness of the mixing angles (see e.g. [7]), they will be neglected. 

Saturating Eq. (0) with the 77 meson state, for the phenomenological part 
of the correlator, we get: 

n = (v(p) \jj\nip + q)){v(p + g) Ijglg) ,^ 

^ (jp + q) 2 - m 2 

The matrix elements entering Eq. (0J) are defined as: 

(V(p)\tf\v(q + P)) = F?(2p+q)p 

(v(P + q)\fM = -if v {p + q)u (5) 

Putting Eq. (jHJ to (jlj), for the phenomenological part of the correlator, Eq. 
(0), we get: 



ll^P, q) = -./„F„ (p + g)2 _ m , (6) 



From Eq. (jUJ), it is seen that there are four different structures that can 
be used to extract the form factor. In this work we will concentrate only on 
the Pfj,p u and q^q v structures. Hence one can write 

U^p, q) = U^q 2 , {p + q) 2 ) Pll q u + U 2 {q 2 , (p + q) 2 ) qil q v + ... (7) 

where . . . stand for the remaining structures and 

Mq 2 ,{p + q?) = -2ifv - ( 5| 2 Q2) 2 

{p + g) 2 - 

Mq\(p + qf) = -if V ( 2 ( g ) 
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where Q 2 = —q 2 . For the invariant amplitudes II i (g 2 , (p + q) 2 ), one can write 
a general dispersion relation in {p + q) 2 in the form: 

/pAs q 2 ) 
ds l — L — + subtraction terms (9) 
s — (p + q) 1 

where the spectral functions corresponding to Eq. (jSJ) are 

p 1 (s,q 2 )=2p 2 (s,q 2 ) = 2zF^(Q 2 )5(s-m 2 TI ) (10) 

Our next task is the calculation of the correlator from the QCD side. If 
the {p + q) 2 and q 2 are large and space like, the OPE in Eq. can be 
expanded near the light cone in terms of the r\ meson wave functions, i.e. the 
matrix elements of the nonlocal operators between vacuum and the r\ meson, 
like (0|grj(/|0) where Tj is an arbitrary Dirac matrix. As we already noted 
that we will take into account twist 3 and twist 4 wave functions beyond the 
leading twist and keep terms up to linear order in the strange quark mass. 
Up to twist 4, the n meson wave functions are defined in the following way 

m 



(V(p)\q(xh^ 5 q(0)\0) = -if vPfi [* due^ x 
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(ri(p)\q(x)i%q(0)\Q) = UN f due^ x VP {u) 



fa(p)|?(aO<V7s?(0)|0) = 

L _ j£\ [paX ^ _ ^ r 1 duel u PX ^ {u) 

" V m r)J-q J ^° 

(r]( P )\q(x)cr a p-f 5 G flu (ux)q(0)\0) = 
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fv m v 
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2 Pm 



T>ai(p±(ai)e 
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if^w 2 



Pp 5V 



px j 



p a g^ 



X(3PfM 

px 



^a^±(ai)e ipz{Q1+UQ3) 
(11) 



ifr,m^(p a x p -p p x a ) J Va^a^e^^ 
After contracting quark fields, the correlator function becomes: 
IV = ]T e i Qy'd 4 xe^(77(p)|%(x)5 i (x)g i (0) + ^(0)^(-x)g i (x)|0) (12) 

i=u,d,s 

where Si(x) is the full propagator of the light quark, and c u ^ = 1 and c s = — 2. 
During the analysis, the following full light quark propagator is used: 



S(x) 



i $ 

2tt 2 x 4 ~ An 2 x 2 



-ig s I 

Jo 



du 



;G ,, v {ux)(J 



16tt 2 x 2 ^ v ' 



ux ll G^ u {ux)Y 



An 2 x 2 



—i 
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where A is the energy cut off separating perturbative and non perturbative 
regimes. 

In order to suppress the contributions of higher states and continuum and 
also to eliminate the subtraction terms in the dispersion relation, the result 
is Borel transformed with respect to — (p+ q) 2 After lengthy calculation for 
the theoretical part of the correlator function, we get the following result 

V6n™(Q 2 ,M 2 ) = 



2{e d + e u - 2e s )^^ J du j dk—^ [B(u) + B{u)\ + q») {kp v + q v ) 

f 

Jo 



«(») 

r -1 e m 2 
+ / du- 



2m s e s 



m 2 / 2 ; M 2 



— [^(w) - v5 ct (m)] (p„qn-p»q v ) 



+2m s e s f ri fj, r] u [(fip(u) - (pp(u)] g^ v 

-^{e d + e u - 2e s )f v [(p v (u) - (p v {u)] 

x (Au 2 p^p u + 2u(p v q lx + p^q u ) - (Q 2 + m 2 u 2 )g^ 

I f ?Ti 2 

+ 8 (ed + e ^ 2es W [A(u) " i(M)] 



X 



(4u 2 PnP u + 2u{p v q il + p M g„) + (Q 2 - M 2 w + m 2 ^ 2 )^) } 



+ dv Va 



e 



'(«! + va 3 ) 2 

x [(« - v){a x + va 3 )p^ t p IM + vp^q u - vp v q^ 

1 /r,m 2 1 

--(e d + e u ~ 2e s ) — j—^— ) 



f Tfl 2 

4(e d + e u - 2e s )^/ [Al^) + Ai(ai)] 



(w - v) (A\\{oii) + ^||(aO) + (V|| (cn) + V||(a, 
4(ai + va 3 ) 2 p^p^ + 2(ai + va 3 )(p u q^ + p^ 
(Q 2 + (ai + ua 3 ) 2 m 2 - («i + t>a 3 )M 2 ) ^ 

\e s mj-^- [vT(a>i) - vT{ai)} (4p M p„ - m 2 ^ 



M 2 

+4(e d + e u - 2e s 



' M 2 



[V±(ai) - V±(aj)] [(ai + va 3 )p IM p I/ + vp^q u + vp v q^ 



(14) 



where M is the Borel parameter and 



s(fc) = - (m 2 /c + Q 



k 

where Q 2 = —q 2 . The contributions of the higher states and continuum are 
subtracted using the quark hadron duality (The details can be found in j^j). 
Hence for the invariant functions in Eq. (JHJ), we get the following results: 



y/6Ut(Q 2 ,M 2 



1 due ~^ J * ( ed + 6u _ 2e s )^ (A(u) - A(u)) 



+2(e d + e u - 2e s 



1 M 2 



u{B{u) + B{u)) 



-2(e d + e u - 2e s )f v [<p n {u) - <Pi,(u))} 



+ dv VaiC 



W 2 



-2(e u + e d - 2e s 



M 2 ai + va 3 
(v - v) (A\\(cti) + A\\(cti) + 2A ± (ai) + 2A±(a 
-V||(oi) + V||(oi) - 2V±(oi) + 2V ± (a 



+lQm s e s 



M 2 («! + TO 3 ) 2 
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{vT(a>i) - vT(a.i)) 



(15) 



P 9 -, s ( fc ) 

V6n 2 ( g 2 ,M 2 ) = 2(e d + e u -2e s ) J -j^ J g duJ o dk^-[B(u) + B(u)) (16) 

3 Numerical Analysis 

In this section, we will carry out the numerical analysis our results. We 
used the following values of the constants appearing in Eqs. (fTHj) and (jTHj) 
0: m a (l GeV) = 150 MeV, f v (l GeV) = 0.156 GeV, and ^(1 GeF) = 

0. 034 GeF 2 . 

In numerical analysis, we systematically take into account all 0(m s ) ~ 
0(m 2 ) effects, which have following sources: a) from the correlator when 
we take p 2 = m 2 ; b) from s-quark propagator, where terms linear in m s 
are retained; c) SU(3) flavor symmetry breaking effects which can appear 
in three different ways: i) via SU(3) violation in normalization factors, i.e. 
fv fn] h) Meson mass corrections to the twist four wave functions; iii) 
via an asymmetry of the wave function with respect to the interchange of 
quark and antiquark. At leading twist two wave function, this leads to the 
appearance of non vanishing odd coefficients of the Gegenbauer expansion, 

1. e. 



<p v (u, jtt) = Quu 



l + 5>n(^/ 2 (2«-l) 



n=l 



where C^ 2 (2u — 1) are Gegenbauer polynomials. In our calculation, we only 
considered ai(/x) and ci2(^) 

From Eq. (|15jh it follows that leading twist two wave function appear 
only in the form ip(u) — ip(u), i.e. they contribute as a result of pure SU(3) 
violating effects. In the numerical calculation, we kept the terms oc ai(//) 
only. The first coefficient was calculated in for K and K* mesons. 

In the present work we have estimated ai(yu) for the rj meson and found that 
a\(fi) = —0.20 ± 0.02. The 0(a s ) corrections to the leading twist two wave 
function has been calculated in [3] in the chiral limit and this correction is 
also included in our calculations. 

For the analysis of the sum rules, one should determine a value of the 
continuum threshold and a region of the Borel parameter M 2 in which the 
predictions do not depend on the value of M 2 . From the analysis of the 
mass sum rules for rj, the continuum threshold so is found to be around 
s — 2.0 GeV 2 . In our analysis we will take s = 1.5 GeV 2 and s = 
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2.5 GeV 2 to study the dependence of our predictions on s . In choosing the 
suitable region for the Borel parameter M 2 , one should keep in mind the M 2 
should be big enough to suppress the contributions of higher dimensional 
operators in the operator product expansion of the correlator function, and 
at the same time small enough to have a considerable suppression of the 
contributions of the higher states and the continuum. In our analysis, we 
will take 0.7 < M 2 < 1.5 in which region the contributions of both the 
higher dimensional operators and the contributions of the higher states and 
continuum are less then 30%. In Figs. (Q) and (J2J), we plot the dependence 
of Q 2 F V (Q 2 ) on the Borel parameter M 2 for Q 2 = 2, 7, 10 GeV 2 and for 
the values of the continuum threshold Sq = 1.5, 2.5 GeV 2 . From the figures 
one sees that the predictions from the p^p u structure are less stable for larger 
values of Q 2 , which is reasonable as the sum rules is expected to break for 
very large values of Q 2 . For smaller values of Q 2 , the predictions are almost 
independent of the exact value of both so and M 2 . 

In Figs. (J3J) and (jH), the dependence of Q 2 F V (Q 2 ) on Q 2 is depicted at 
the value of the Borel parameter M 2 = 1 GeV and for two different values of 
the continuum threshold Sq = 1.5 and 2.5 GeV. As can be seen also in these 
figures, the predictions obtained from the structure p^p u is more robust to 
variations in the continuum threshold sq then the predictions obtained from 
the q^q u structure. A fit to the form factor F V (Q 2 ) can be written in the 



From the analysis of the sum rules obtained from the p^p v structures, for 
the constants we obtained the following values: ^(0) = 0.17 ± 0.02, 
6i = (8 ± 2)10~ 2 , and b 2 = (4 ± 1)10~ 3 . 
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Figure 1: The dependence of Q 2 F V (Q 2 ) on the Borel parameter at Q 2 
2, 7, 10 GeV 2 for the (a) p^p u and (b) structures 
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Figure 2: The same as Fig. (0) but at sq = 2.5 
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Figure 3: The dependence of Q 2 F rj (Q 2 ) on the Q 2 obtained from the structure 
V^Vv at M 2 = 1 GeV 2 and at s = 1.5, 2.5 GeV 2 
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